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Abstract — In this letter, we investigate the resource allocation 
for downlink multi-cell coordinated OFDMA wireless networks, 
in which power allocation and subcarrier scheduling are jointly 
optimized. Aiming at maximizing the weighted sum of the 
minimal user rates (WSMR) of coordinated cells under individual 
power constraints at each base station, an effective distributed 
resource allocation algorithm using a modified decomposition 
method is proposed, which is suitable by practical implemen- 
tation due to its low complexity and fast convergence speed. 
Simulation results demonstrate that the proposed decentralized 
algorithm provides substantial throughput gains with lower 
computational cost compared to existing schemes. 

Index Terms — Resource allocation, multi-cell network, opti- 
mization theory 

I. Introduction 

OFDMA has been widely accepted as a promising multiple 
access technique for future mobile communication systems. 
Recently, OFDMA-based multi-cell coordination becomes a 
hot research topic, which could provide superior performance 
over the traditional single-cell processing network through 
joint signal processing among the involved based stations 
H). Moreover, the performance of the multi-cell coordination 
network can be further improved by employing multi-cell 
resource allocation. 

However, the traditional centralized multi-cell resource allo- 
cation approach is too complicated to be practical because of 
the large number of parameters and constraints Q. Therefore, 
various distributed algorithms have been proposed to reduce 
the computation complexity and increase the scalability, such 
as the distributed algorithm for the resource allocation using 
game theory in and the one based on the Lagrange duality 
method in 0). 

In this paper, a novel distributed resource allocation algo- 
rithm is proposed for multi-cell OFDMA networks. Different 
from existing works in which weighted sum of user rates 
(WSR) are usually used, the weighted sum of the minimal user 
rates (WSMR) of coordinated cells is taken as the performance 
criterion, which should be maximized subject to individual 
power and subcarrier allocation constraint at each BS. An 
iterative algorithm is proposed to solve the problem, which 
optimizes the subcarrier scheduling and the power allocation 
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alternatively. For the power allocation subproblem, motivated 
by a natural decomposition of the optimality conditions of the 
original problem, we propose a novel distributed algorithm. 
Unlike traditional Lagrangian based algorithms, the proposed 
technique does not need to solve subproblems, which results 
in computational savings and fast convergence. Furthermore, 
in the proposed procedure the central agent only distributes 
information and checks the convergence condition without the 
need to update information as J4), which makes the scheme 
simpler. 

II. System Model And Problem Formulation 
A. System Model 

In this paper, we investigate the distributed resource alloca- 
tion in the downlink of a cellular OFDMA network consisting 
of M cooperative cells. In the m th cell, BS m serves K m 
users and the active user set is denoted as U m . The channel 
power coefficient of subcarrier n from BS / to user u in the 
m th cell is denoted by g l u m n . If I ^ m, g u m n is related to 
the interfering channel power from BS I , and g l umn denotes 
the desired channel power coefficient from BS m if I = m. In 
addition, a\ m n is the power of the additive white Gaussian 
noise at user u in cell m on subcarrier n. The number of 
OFDMA subcarriers is N and each subcarrier is allocated to 
only one user exclusively in each cell. It is also assumed that 
all channel state information (CSI) is perfectly known at each 
BS. 

The maximum transmit power of BS m is denoted as 
Pm,max, and P m ,n represents the power allocated to subcar- 
rier n by BS m. For convenience, {P m ,n}m=i is stacked 
into a M x 1 vector p„ = [Pi, n , . . . , Pm,e] T . an d then 
{PnJ^Li is stacked into a M X N matrix P referring to 
power allocation in following sections. As for the subcarrier 
allocation notations, we define a binary variable A u>m n , which 
indicates that subcarrier n is allocated to user u in cell m 
if Ai,m,« = 1. Stack {A Utm . n }% =1 into a JV x 1 vector 
A u , m = [A u . mA , . . . , A u , m . N ] T , and merge A u<m of all users 
in all cells to a W x K m x M subcarrier allocation matrix A 
which indicates how subcarriers are assigned among all users. 

Then the information rate of user u in cell m can be 
expressed as a function of A u>m and P 

Pu,m ( A u?m j P) — y A uran ■ -Rii,m,n(Pn) 
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where R u . m ,n is the achievable rate of user u on subcarrier n 
in cell m, and can be calculated as 



Ru,m,n(Pn) — In [ 1 



P, a m 

1 m,ny u ,m,n 



(2) 

in the unit of nats per OFDM symbol, where 7 Ujm .„(p„) is the 
signal to interference plus noise ratio (SINR), and V represents 
the signal to noise ratio gap between the adopted modulation 
and coding scheme and the one achieving capacity. 

B. Problem Formulation 

For resource allocation, sum capacity is the most widely 
used performance criterion. However, since OFDMA multi- 
cell networks with multiple users in each cell, maximizing the 
sum capacity of a multi-cell network may result in a serious 
performance imbalance. It is because that the users and cells 
with better channel conditions will be allocated with much 
more resources and those experiencing worse channels may be 
sacrificed. In order to overcome this kind of problems fairness 
among the users and cells must be taken into account. 

Therefore, an objective function called weighted sum mini- 
mal user rate (WSMR) is used. As shown in 0, maximizing 
the minimal user rate in a single cell could provide the 
maximum fairness among the users in this cell. Therefore, 
in cooperative multi-cell systems, by introducing different 
priorities to the involved cells, WSMR can be employed 
to provide fairness to multiple users as well as cells. The 
objective function in this problem is given by 
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ug(7 m 



where uj m > represents the weight assigned to cell to's 
minimal user rate. In particular, increasing oj m leads to a 
higher resource allocation priority assigned to the users in cell 
to. Based on (0), the resource allocation problem is formulated 

as 



max f(A.P) 
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Note that the minimization operation in (HJi prohibits the 
objective function from being differentiable. Thus, an auxiliary 
optimization variable R m is introduced, and the optimization 
problem is reformulated as follows 



max 
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Obviously, the optimization problem © is a constrained 
nonlinear optimization program with both integer and continu- 
ous variables. Furthermore, the problem is not convex and thus 
generally speaking it is difficult to directly solve the problem. 
Following a similar logic as those in O, Q, an effective algo- 
rithm is proposed in the following, which optimizes subcarrier 
scheduling and power allocation alternatively. 

III. The Proposed Algorithm 

A. Proposed Decomposition Power Allocation Algorithm 

Without loss of generality, the power allocation optimization 
is carried out first at each iteration. The subcarrier allocation 
variables A can be directly removed from (O as they are 
assumed to be already computed. For notational simplicity, all 
superscript of variables are omitted in following description. 
Thus, the power allocation optimization subproblem can be 
expressed as 



max 
s.t 
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Based on the definition of R u ,m.n(Pn) in it is obvious 
that (O is a nonconvex optimization problem. For multi-cell 
OFDMA networks, the optimization problem ^ usually has 
high dimensions. Although centralized optimization methods 
can be used to solve the problem, they are impractical due 
to the high complexity. Therefore, instead of centralized opti- 
mization algorithms, we focus on distributed power allocation 
schemes using decomposition techniques. 

Various decomposition algorithms have been developed. 
The frequently used decomposition algorithm in wireless 
communications is Lagrangian relaxation procedure 0, (9), 
and its variant relaxation techniques based on augmented La- 
grangian functions iflOl . 1111 . However, Lagrangian procedures 
may present drawbacks in some cases, such as difficulties 
to converge to an optimal solution for the global system 
(in the absence of convexity assumptions), uncontrollable 
convergence rates that depend on the correct choice of the 
values for several parameters which are difficult to update, 
and the requirement of the intervention of a central agent to 
update this complicated information. 

To overcome these drawbacks, a novel decomposition algo- 
rithm which is based on the decomposition of the optimality 
conditions of the global problem © is proposed in this paper 
0. The proposed decomposition algorithm improves both the 
computational efficiency and implementation simplification 
compared to previously mentioned algorithms. In the follow- 
ing, the decomposition methodology is presented in detail. 

First of all, the optimization problem §6^ can be written in 
a compact form as follows for convenience 



EM 
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s.t h(x-L,...,x M ) <0 g m (x m )<0 
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where x m is a vector for cell m, that contains power allocation 
vector p m and the minimum user rates R m within cell to. 

The first constraint in (0 is known as complicating con- 
straint, which represents the rate constraint in (6). These 
equations contain variables and parameters from different cells 
and prevent each system from operating independently to each 
other. On the contrary, the second constraint is the power 
constraint for each BS and thus only related to one single 
cell. Therefore, in order to decompose the global optimiza- 
tion problem (0 into local subproblems for each cell, those 
equations of rate constraints are removed from 0. Thus, the 
problem is equivalent to 



is formulated as 



,/m(aVi)+VV \^ n h m {x-i,...,x M ) 
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s.t h m (x 1 ,...,x M ) < g m {x m ) < 
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where the first constraint has been separated into different 
cells compared with that in 0. The dual variable vector 
corresponding to the first constraint is denoted by A m . 

Fixing the values of all variables and multipliers (indicated 
by over line) except those in cell m, (0 reduces to 

A, hi(x m ) 

S.t X = [X\) . . . , 33 m _i , , X m _|_i , . ■ • , 3?m] 

M* m ) < g m (x m ) < (9) 

where k — J2t=i i=£ m //(*/) i s a constant. Similarly, the 
reduced problem ^ can be reduced for every cell (to = 
1,...,M). 

To prove that the proposed decomposition method is based 
on the solutions of these reduced cell subproblems, we revisit 
the first-order optimality conditions of the problem which 
is written as 

EM T t * * 

. ^x m hm{xl, . . . , X* M )\ m 
m—1 

+ Vl m g m {x* m )n* m = 

h m (x*, . . . , x* M ) < 0, h m (x m ) T \ m = 

Kn > °> 9 m «) < 0, 9 m {x* m ) T tf n = 
H* m >0 m=l,...,M (10) 

These conditions have been constructed using the optimal 
values x* m , X* m and fi* n which are assumed to be known. The 
values A* n and the values /j, m are the optimal Lagrange multi- 
pliers associated with the two constraints in 0, respectively. 
If the first-order optimality conditions of every cell reduced 
subproblem (0 (m = 1, . . . , M) are stuck together, it can be 
observed that they are identical to the first-order optimality 
conditions ( TTOb of the global problem (0. Thus, it is obvious 
that the decomposition mythology holds. 

Based on the previously proposed decomposition method- 
ology, removing constants and substituting the objection func- 
tion and constraints of (0 into (0, the original power alloca- 
tion problem can be decomposed into M subproblems, which 
are computed in each cell independently. The to" 1 subproblem 



max uj m R m + y^ X u ,i(R, 
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where p m = (pi, • • • , p m _i, p m , Pm+i, • • • , Pa/), and A„ >m 
is the optimal Lagrange multiplier which guarantees that user 
rate is larger than minimal value R m in each cell. 

A summary of the proposed decomposition algorithm for 
power allocation optimization is as follows: 

Algorithm 1: Power Allocation Optimization Algorithm 
Initialize: Initialize the iteration counter t = 0, and each 
cell (m = 1, . . . ,M) initializes its variables 
Pm,R m and parameters A UiTO . 

repeat 

1. Each cell carries out one iteration for its 
corresponding subproblem QT}, and obtains 
search directions Ap m , AR m , AA„, m ; 

2. Each cell updates its variables and parameters 

Pm Pro + Ap m , R m R rn + AR rn , 
Au,m ^ X u ^ m -\~ AA Ujm? 

3. t = t + 1; 

until || vcc(P* - P*" 1 ) || < * , or t = T; 

Return P* 

The search directions, Ap m , AR m , AX u , m , for subproblem 
(fTTT > can be computed independently of each other, allowing 
a parallel and distributed implementation. This step requires a 
central agent to coordinate the process, which receives certain 
information (p m , R m , X u . m after each iteration) from all cells 
and returns it to the appropriate cells. It can be noted that the 
information exchanged between the areas and the central agent 
is little. Unlike other decomposition algorithm, in the proposed 
algorithm the central agent only distributes information and 
checks the convergence condition. It does not need to update 
any information, because this information is updated by the 
areas of the system, implying a simpler process. 

The main difference between the Lagrangian relaxation 
algorithm and the proposed decomposition one is that La- 
grangian relaxation adds all the complicating constraints into 
the objective function. Therefore it needs auxiliary procedures 
to update the Lagrange multipliers. On the contrary, the 
proposed technique does not need any procedure to update 
the multipliers because this updating is automatic and results 
directly from the foreign optimization problem. For example, 
the Lagrangian multipliers A; (I ^ to) included in the 
objective of cell to are obtained from cell I 's optimization 
by keeping its own complicating constraints. What's more, the 
proposed approach has the advantage that convergence prop- 
erties do not require an optimal solution of the subproblems 
at each iteration of the algorithm. It is enough to perform 
a single iteration for each subproblem, and then to update 
variable values. As a consequence, computation times can be 
significantly reduced with respect to other methods that require 
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the computation of the optimum for the subproblems in order 
to attain convergence. 

B. Subcarrier Allocation Optimization 

When the power allocation is fixed, the remaining optimiza- 
tion problem is to find the optimal subcarrier allocation. It can 
be directly decomposed into M subproblem, each of which 
corresponds to problem of maximizing the minimal user rate 
in cell to. The to*' 1 is a Mixed Integer Linear Problem (MILP) 
problem involving both integer and continuous variables, it 
can be solved by various algorithms such as exhaustive search, 
implicit enumeration method and branch-and-bound algorithm 

m. 

IV. Simulation Results 

In this section, the performance of the proposed distributed 
resource allocation algorithm is investigated. The downlink of 
a cellular OFDMA system is considered with three coordinated 
cells and 32 subcarriers. The radius of each cell is 40m 
and there are two users randomly located in each cell. To 
focus on the performance of proposed algorithm, all weighting 
coefficients assigned to each cell are assumed to be identical 
and keep constant. In addition, we set T = OdB, a 2 ^ m n = —60 
dBw and *i = * 2 = *3 = 10 _1 . 500 independent channel 
realizations are carried out to obtain the final results. 

FigOJ compares the average WSMR performance of the 
proposed algorithm and the one using Lagrangian relaxation 
method (LR) as a benchmark power allocation algorithm. The 
average WSMR at initialization which use uniform power 
allocation and even subcarrier allocation is also shown as a 
reference. It can be seen that the performance of the proposed 
distributed resource allocation algorithm is better than that 
of the algorithm using traditional Lagrangian decomposition 
method. Figf2] demonstrates the convergence performance of 
different algorithm for power allocation. Unlike the slow and 
oscillating behavior of LR procedure, proposed algorithm 
converges rapidly and stably. It can be concluded that the 
proposed decomposition algorithm are more effective than the 
traditional Lagrangian based algorithms. 

V. Conclusions 

In this letter, we proposed a distributed algorithm for joint 
resource allocation in a coordinated multi-cell OFDMA net- 
work. The performance criterion termed as WSMR has been 
introduced to guarantee the fairness among multiple users and 
cells. Targeting at maximizing WSMR, an iterative algorithm 
was proposed to solve the joint optimization of power allo- 
cation and subcarrier scheduling. At each iteration, the power 
allocation is updated by applying a modified decomposition 
methodology which has low complexity and fast convergence 
speed. On the other hand, the subcarrier allocation is updated 
by solving a MILP. It is shown that the proposed distributed 
resource allocation algorithm achieve a better WSMR perfor- 
mance compared compared with the traditional decomposition 
algorithms. At the same time the proposed power allocation 
scheme also provides a better convergence performance than 
the Lagrangian algorithm. 
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Fig. 1. Average WSMR of proposed algorithm and the algorithm whose 
power allocation use Lagrangian relaxation method. 



^ f ■ 1 

<S> A - i 



— •— Proposed algorithm 
■ -»- ■ Lagrangian Relaxation 



Iteration Index 

Fig. 2. Convergence performances of proposed decomposition algorithm and 
the Lagrangian relaxation method for power allocation. 



